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Abstract. We investigate the sparsity of the Gabor-matrix representation of Fourier 
integral operators with a phase having quadratic growth. It is known that such 
an infinite matrix is sparse and well organized, being in fact concentrated along 
the graph of the corresponding canonical transformation. Here we show that, if 
the phase and symbol have a regularity of Gevrey type of order s > 1 or analytic 
(s = 1), the above decay is in fact sub-exponential or exponential, respectively. We 
also show by a counterexample that ultra-analytic regularity (s < 1) does not give 
super-exponential decay. This is in sharp contrast to the more favorable case of 
pseudodifferential operators, or even (generalized) metaplectic operators, which are 
treated as well. 



1. Introduction 
We consider Fourier integral operators (FIOs) in the reduced form 

(1) Tf{x)= [ e 2 ™*( x rta(x,r,)f(ri)dr, 



of the type of those in [TJ, namely the amplitude a(z), z = (x, rf), is in Sq , i.e. d"a(z) 
is bounded for every a and the real- valued phase function &(z), satisfying the standard 

(2) 

nondegeneracy condition, belongs to Sqq, i.e. d^{z) is bounded for |a| > 2. Such 
FIOs represent the propagators at a fixed time t > 0, for the Schrodinger equations 

(2) D t u + a w (t,x, D x )u = 0, u\ t=0 = f(x), 

(2) 

with real- valued Hamiltonian a(t, x, £) belonging to Sq q uniformly in t, see for example 
Tataru [37] and Bony [HE]. 

In [12\ [37] it was proved that the Gabor matrix representation of T is concentrated 
along the graph of the canonical transformation x determined by and provides 
optimal sparsity. The Gabor representation was then used to discuss the bounded- 
ness properties of T, cf. [9l [TT] , and define Wiener algebras of global FIOs containing 
operators of type ([1]), see [TO] . 

In the present paper, cf. Section 3, we shall present a stronger sparsity result, with 
exponential decay, for the case of analytic- Gevrey functions, namely when we have in 
([TJ for some s > 1, 

(3) \d a <P(z)\ < C^(a\) s , «GN M , |a| >2, z£R 2d 
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and similarly for the amplitude: 

(4) \d a a(z)\<C M (a\) s , a G N M , z G R 2d . 

As a side result, we shall deduce boundedness of T on Sg(Mr), Gelfand-Shilov spaces 
(basic definitions and properties for these spaces are recalled in the preliminary Section 
2). 

We shall not give explicit applications to the general Schrodinger equation ([2]) in 
the present paper. As a matter of fact, a precise version of [lj in analytic-Gevrey case, 
i.e. when the estimates of the type ([3]) are satisfied by the Hamiltonian a(t,x, £), is 
missing in the literature as far as we know. Note however that there is a number of 
papers where the Schrodinger propagators are treated in the analytic framework under 
decay assumptions for a, <£, cr, see for example [2U Ell EHl EHJ EU [32]; cf. also [BJ [7] 
and [30, Chapter 6] for standing wave solutions. 

Let us also mention the reach literature concerning the different case of the Hdrman- 
der's FIOs [25] ; i.e. positive homogeneity of degree 1 with respect to r/ for $(x,r/) and 
corresponding decay estimates in ([3]), ([!]), mainly addressed to the study of the hyper- 
bolic equations. For such FIOs in the analytic-Gevrey category see the bibliography 
of [33] concerning the intensive production of the years '80-'90. The researches in this 
area are indeed extremely active also nowadays, with applications to weakly hyperbolic 
problems in Gevrey classes. 

In the above mentioned literature, the analytic regularity s = 1 is regarded as opti- 
mal result. Instead, when dealing with the Gelfand-Shilov classes Sg(M, d ), it is natural 
to question whether we can go beyond the barrier s = 1, getting super-exponential 
sparsity and ultra-analytic regularity, i.e. boundedness in Gelfand-Shilov spaces for 
1/2 < s < 1. As we shall clarify in Section 4, this is possible if and only if the phase 
function &(x,r/) is quadratic in x,rj. Such propagators are obtained from ([2]) when 
a(t, has quadratic principal part in (£,£)• The corresponding operators T in ([T]), 
with amplitude satisfying (JI]) for 1/2 < s < 1, are studied in Section 5. 

In the second part of this introduction we want to give a short presentation to 
Gabor frames, addressing to non-expert readers. Generally speaking: paradigm of 
the applications of harmonic analysis to the study of operators and function spaces is 
the decomposition/reconstruction into "wave packets": Fourier series, wavelets, para- 
products, etc., see the survey work [40] and also [35] for applications to dispersive 
equations and the restriction theorem. We may say that every class of symbols, i.e. 
every class of partial differential equations, requires a corresponding partition of the 
phase space into wave packets, cf. [17] . Consequently, we may represent the propaga- 
tor as an infinite matrix, and the chosen partition works effectively for the problem 
under investigation if the matrix is sparse and well-organized. This means that the 
propagator re-arranges the wave packets with a controlled number of overlapping of 
supports, granting continuity on function spaces. 

For wave equations and Hormander's FIOs let us mention the pioneering work of 
Cordoba and Fefferman [15], the second dyadic decomposition of Seeger, Sogge and 
Stein [33] and the phase space transform of Tataru and Geba [38] ; see also Tataru [36] 
for applications to wave equations with non-smooth coefficients. 
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Numerically stable treatments of Hormander's FIOs were carried on by Candes, 
Demanet [H[5] and Guo, Labate [25j, wave packets being represented by curvelets and 
shearles. 

Gabor frames, used initially for problems in Signal Theory and Time-frequency 
Analysis, cf. |T6j UHl |22] , turn out to be the correct setting for Schrodinger propagators, 
at least when in the Hamiltonian the space variables x and their duals £ play symmetric 
role, as we have in the Sq class. This means that the micolocal propagation of 
singularities is identified by the canonical transformation, modulo errors which we may 
estimate parithetically in the x and £ variables, see [21 El O HOl Ell E21 EZ] , mentioned 
before. 

To be definite, let us recall some basic definition. Let A = AI? d with A G GL(2d,R) 
be a lattice of the time-frequency plane. Consider the time-frequency-shifts 

(5) gx = gx lM = e 2niX2X g(x - AO, A = (A 1; A 2 )gA. 

The set of time-frequency shifts Q(g,A) = {gx : A G A} for a non-zero g G L 2 (M, d ) 
is called a Gabor system. The set Q(g, A) is a Gabor frame, if there exist constants 
A, B > such that 

(6) A\\ff 2 < £ \(f,g x )\ 2 < B\\f\\l V/ G L 2 (R d ). 

AeA 

If © is satisfied, then there exists a dual window 7 G L 2 (M d ), such that ^(7, A) is a 
frame, and every / G L 2 (M. d ) possesses the frame expansions 

/ = ^2(f,g\)i\ = ^2(fn\)9x 

AeA AeA 

with unconditional convergence in L 2 (M. d ). 

The Gabor decomposition of an operator T is then as follows: 

T f{x) = X) X] ( T #A,5/i) c a7m> with c a = (/, 7a}- 
/xeA AeA s Ty ' 

So the action of the operator T above can be read on the coefficient space as 

{ca}a g a ^ { ^(Tgxid^cx}^, 
AeA 

i.e. it is represented as the infinite matrix {M^^j^.AeA = {(Tgx, 5 f /x)}^,AeA ! which we 
call the Gabor matrix of T. 

We can now describe the main results of the paper. Let T be defined as in ([TJ), with 
non-degenerate phase function $ satisfying (|3|) and amplitude a satisfying (j3J). Let 
X '■ R 2d — > R 2d be the canonical transformation associated to 3>. In the generic case 
s > 1, fix a window g G S s s j 2 (M. d ). Then for some e > 

(7) |M M |<exp(-e|/,- X (A)r/ s ); 

see Theorem 13.31 below. Besides, if $ is quadratic then ([7]) keeps valid for s > 1/2, for 
any choice of the window g in Sg(M. d ), see Theorem 15.31 The Gaussian, belonging to 



4 



ELENA CORDERO, FABIO NICOLA AND LUIGI RODINO 



1/2 A 

Sy 2 (M ) would work as window in any case. Sparsity and boundedness follow easily, 
see Propositions 13.41 and 13.51 

A class of counterexamples to the validity of (J7]) when s < 1 and $ is not a quadratic 
polynomial is given in Proposition 14.11 

2. Preliminaries 

2.1. Notations. The Schwartz class is denoted by <S(M <i ), the space of tempered distri- 
butions by S' (R d ). We use the brackets (/, g) to denote the extension to <S'(R ) x <S(K ) 
of the inner product (f,g) = J f(t)g~(Jjdt on L 2 (R d ). 

The Fourier transform is normalized to be f(rj) = J-f{rj) = f f{t)e~ 2mtr) dt. 
Translation and modulation operators, T and M are defined by 

T x f{-) = f{--x) and M x f(-) = e 2mx -f(-), x G M. d . 

The following relations hold 

(8) M y T x = e 27Tixy T x M y , (T x fJ = M- X f, (M x f}=T x f, x, y G R d , /, g G L 2 (R d ). 

Throughout the paper, we shall use the notation A < B to express the inequality 
A < cB for a suitable constant c > 0, and A x B for the equivalence c~ 1 B < A < cB. 

The letter C denotes a positive constant, not necessarily the same at every appear- 
ance. 

2.2. Gelfand-Shilov Spaces. Specially in Applied Mathematics, it is of great interest 
to quantify the decay of functions at infinity, and the Schwartz class reveals to 
be insufficient for this. The so-called Gelfand-Shilov type spaces, introduced in [20] 
turn out to be very useful. Let us recall their definition and main properties; see e.g. 
[201 EQ] for details. 

Definition 2.1. Let there be given s,r > 0. A function f € <S(M d ) is in the Gelfand- 
Shilov type space S^(M. ) if there exist constants A,B>0 such that 

(9) \x a B p f{x)\ < A^B^(a\) r (pi) s , a,/3eN d . 

The space ^(IR^) is nontrivial if and only if r + s > 1, orr + s = l and r, s > 0. So 

1 /2 A 

the smallest nontrivial space with r = s is provided by Sy 2 (R ). Every function of the 
type P(x)e~ a \ x \ 2 , with a > and P(x) polynomial on ~R d , is in the class S^(R d ). We 
observe the trivial inclusions S^IR^) C S^(M. d ) for s± < S2 and r\ < r2- Moreover, if 
/ G Sr(M. d ), also x s d 1 f belongs to the same space for every fixed (5,7. 

The action of the Fourier transform on S'^.(M d ) interchanges the indices s and r, as 
explained in the following theorem. 

Theorem 2.2. For f G S(R d ) we have f £ S°(R d ) if and only iff G S r s (R d ). 

Therefore for s = r the spaces are invariant under the action of the Fourier 

transform. 

We shall also need the following analyticity property of functions in S^.(R d ), when 
s < 1. 
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Theorem 2.3 ([301 Proposition 6.1.8]). Assume f G S£(R d ), < s < 1, r > 0. Then 
f extends to an entire analytic function f(x + iy) in C d , with 

(10) |/(x + *lO| <e- e M* +% 1 T=? ) x G M d , y G M d , 
where e and 5 are suitable positive constants. 

Let us underline the following property, which exhibits two equivalent ways of ex- 
pressing the decay of a continuous function / on R d . This follows immediately from 
[301 Proposition 6.1.5], see also [13, Proposition 2.4], where the mutual dependence 
between the constants e and C below was shown. 

Proposition 2.4 ( |30l Proposition 6.1.5], [131 Proposition 2.4]). Consider r > and 
let h be a continuous function on R d . Then the following conditions are equivalent: 

(i) There exists a constant e > such that 

(11) \h(x)\<e~^, x G R d . 

(ii) There exists a constant C > such that 

(12) \x a h(x)\ < C H (a!) r , x£R d ,aeN d . 

frdY 

Indeed, assuming (jlip , then f)12|) is satisfied with C = ( — j . Viceversa, (|12p implies 

(jlip /or any e < r(dC)~~- . Also, the constant implicit in the notation < in (lllj) 
depends only on the corresponding one in (j!2p and viceversa. 

2.3. Time-frequency characterization of Gelfand-Shilov spaces. Consider a 
distribution / G <S'(R ) and a Schwartz function y G <S(M rf )\{0} (the so-called window). 
The short-time Fourier transform (STFT) of / with respect to g is 

V g f(x, n) = (/, g x , n ) = [ e 2 ^gjt^x)f(t) dt (x, n) G R 2d . 

The short-time Fourier transform is well-defined whenever the bracket (•, •) makes sense 
for dual pairs of function or (ultra-) distribution spaces, in particular for / G ^'(IR^) 
and g G S(R d ), f,g G L 2 (R d ), or / G (S°)'(R d ) and g G S°(R d ) (see [22] for the full 
details). 

The following inversion formula holds for the STFT ([22, Proposition 11.3.2]): as- 
sume g G 5(M d ) \ {0}, / G L 2 (R d ), then 

(13) / = TTm2 / Vgfix^M^gdxdr]. 

\\9h Jm. 2d 

Finally, we have the following characterization of Gelfand-Shilov functions; cf. [H 

[THESES]. 

Theorem 2.5. If s > 1/2, 

(14) f,geS s s (R d ) ^V g f eS s s (R 2d ); 
ifge S s s (R d ), then 

(15) / G S s s (R d ) <=> \V g (f)(z)\ < e^ 1/a , z G R 2d , for some e > 0. 
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3. Exponential sparsity of the Gabor matrix representation 

The Fourier integral operator T with symbol (or amplitude) a and phase on 
is formally defined in ([I]). The phase function $>(x, 77) is smooth on M. 2d , and fulfills 
the estimates 

(16) \d a <S>(z)\ < C^(a\) s , aGN M , \a\ > 2, zeR 2d , 
for some C > 0, s > 1, as well as the nondegeneracy condition 

(17) I det d^(x,rj)\ > 6 > 0, (x,??)elR M . 
The symbol a on R 2d satisfies 

(18) \d a a(z)\ <M{z)C^(al) s , a G N 2d , z G R M , 

for the same s as in (|16j) and some C > 0, and some continuous weight M > in 
We assume here that M is temperate, in the sense that 

(19) M(z + w) < (z) N M(w), z,w£R 2d , 
for some N > 0. 

We also denote by x '■ ^ M — * ^ 2<i the canonical transformation defined by $, i.e. 

J y = V r ,<f>(x,ri) 
\ £ = V x $(x,ri). 

The canonical transformation % enjoys the following properties: 



(20) (X,Z)=X(V,V) 



(i) x '■ ^- 2d ~ * ^- 2d is smooth, invertible, and preserves the symplectic form in M. 2d , i.e., 
dx A d£ = dy A dr\\ x is a symplectomorphism. 

(ii) For 2 = (y 5 7?), 

(21) |9? X (*)|<cW(a!)', |a|>l; 
(mj There exists 5 > such that, for = xiUiV)) 

(22) |det— (y,^^. 
We need a preliminary lemma. 

Lemma 3.1. Let s > 1 and <£>(z) a reaZ smooth function in M. d satisfying the estimates 

\d a y{z)\ < C^ +1 (a\) s (z} 2 , a G N d , G R d , 
/or some constant C > 0. ITien /or t/ie same constant C it turns out 



|9 a e^| < (d»- 1 2 d+1 C7 2 )H ^ (^)%) 2i , |a| > 1, z G 
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Proof. By the Faa di Bruno formula (see, e.g., [El page 16]) and the hypothesis we 
have, for |a| > 1, 

M 1 , 

|fl"e w l < E - E r 1 , l g7 V(^)l • • • l^(*)l 

j = l J 7i+...+7j=a ' J 



l7fcl>l 

< E ^ E - r ^c^i+-i^w( 7l ! . . . 7j ir(z) 

7 ! 71 ! . . . 7, ! 

j = l J 7 1 +...+ 7j =a fi >J 

l7fel>l 

= E^w|(^- £ (Tii.-.T.ir 1 . 

j'=l 7l+---+7j=<* 



2j 



l7fel>l 

Now will verify that 

Id! 

(23) 7i + ... + 7j = a, |7fe| > 1 71! . . . 7,! < — -. 

This then gives the desired conclusion, taking into account that s > 1, |a|! < d' a 'a! 
and 

e 1 - n + 3 ~ x ) - 2 |a|+d( ^ i) < 2( d+i )i a i. 

7 1 +...+7 J -=a fc=l \ 3 J 

l7fcl>l 

It remains to prove (|23p . We argue by induction on j. If j = 1 it is obviously true. 
Let therefore j > 2 and assume that (j23[) holds for j — 1 factors. Then 

71-... 7,-1-7,.- (J _ 1}! 7,-_ 17,1- 

_ 3 (H ~ l7jl)!|7j|! 

j! ' H ' (H-l)! 

Since j < |a| the derided estimate in (|23p therefore follows if we prove that the last 
fraction is < 1. But this is clear because 

(H-| 7j |)!| 7j |! _ | 7j | | 7j -|-i l 7j -|-2 ^ 2 

(|a| — 1)! \a\ — 1 \a\ — 2 |a| — 3 |ck[ — |7j| + l 

and in this product each fraction is < 1: indeed, j > 2 and |ti| > 1 imply |7j| < |a| — 1 
and therefore \ — k < \a\ — 1 — k, for < k < |7j| — 2. □ 

Remark 3.2. Let us observe that the Faa di Bruno formula, combined with the for- 
mula (|23l) . gives a cheap proof that Gevrey classes are stable by functional composition, 
with precise estimates for the constants; we omit the details. 

s > 1, 

above. Assume g € S S J^(M^) Then there exists e > such that 



Theorem 3.3. Let s > 1, and suppose the phase <3? and symbol a satisfy (|16p -(|T9 
above. Assume g € S^^iW 1 ) Then there exists e > such that 

(24) \(Tg u ,g v }\ < M(vi,u 2 )eyip (-e\v-x{u)\ 1/s ), u = (u 1 ,u 2 ), u = («i,« 2 )G 



2<Z 
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Proof. A direct computation based on (|8|) (see e.g. the proof of [JTJ Theorem 3.1]) 
shows that 

(Tg u ,9v) 

e a 9 rf[*(«+« 1 ,fl+« a )-(« 3 ,«0.(*+wi I ij)] o .( a . + VuV + U2 )g(x)g(r } ) dxdr] 
By performing a Taylor expansion of around (vi,u 2 ) we obtain 

(25) \(Tg u ,g v }\= [ e ^.%wH*.«i))' e ^*2,(«i,« l )W ( r( 2 + m 2 ))G(^) cfe 
where Cr(z) = G(x,rj) = g(x) ®'g(rj), and 

(26) ^,(„ 1) « 2 )W = 2V / (l-t^t^i.^l + k)/, 2 = (z,7?), 

| Ct |=2- / ° ° ! 

is the second order remainder in the Taylor formula for <& at {vi,u 2 ). Observe that 
(|16p implies the estimates 

(27) \d^ 2i{viiU2) (z)\ < C^a!)^} 2 , a G N M s € M M 

uniformly with respect to V\,u 2 G W 1 . 

Now, it is proved in [IT] , Lemma 3.1] that 

[V*$(vi,lt2) - («2,«l)| > b - X(u)\. 

Hence, it is sufficient to prove that 

\(Tg u ,g v )\ < M(ui,«2)exp ( - e|V«*(ui,«2) - (^,«l)| 1/s ), u,v G M 2d 

Using the formula ([25]) for the left-hand side, we are reduced to proving that the 
function 

h VliU2 (uj) := f e 2muJz e 2m ^i^ iz) a{z + (v 1 ,u 2 ))G{z)dz, uj G M. 2d 
satisfies the estimates 

\ h v u u 2 (u)\ < M(v 1 ,u 2 )ex 1 p ( - e\u\ 1/s ) 
or equivalently, by Proposition 12.41 

(28) \Lj a h VltU2 (u)\ <M(v u u 2 )C^(a\y, a G N 2d , u G R 2d . 
Now, repeated integrations by parts and Leibniz formula give 

(29) \0J a h Vl , U2 (u)\ < (27r)-l a l / (z)-*^ 1 ***" £ 

x ( z ) 2d+1 d^e 2ni ^i^)( z )d^a{z + {v u u 2 ))d^G{z)dz 
Let us estimate the three derivatives above. By (|27p and Lemma 13. II we have 

1/3x1 l 

(30) |^e 27r ^,(, 1 ,„ 2) W| < C \fh\ ^y(z} 2 ^, 1/3x1 > 1. 
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Using ([15]) and (fT9[) , the derivatives of the symbol can be controlled by 

\d^a(z + {v l ,u 2 ))\ < M{v x ,u 2 )C\^\{^ a {z) N . 

Hence, for > 1, 

(31) \{z) 2d+l d^e 2m ^^ z )d^a{z + {v u u 2 ))d^G{z)\ 

m ,R Ws 

<M(v u n 2 )C^ + ^K^r(z) N+ ^J2{^r) (z) 2 '\d^G(z)\. 

3=1 J ' 

Now, by Theorem E2 we have g € S s £(R d ), so that G = g®g € S s / / l(M. 2d ) . This 
gives 

^+2^1+2^/^)1 < C N+2d+i+2j+m^ N + 2 d+l + 2i)!) s / 2 (/3 3 !) s/2 
where we used the formula (m + n)\ < 2 m+n m\n\ and Stirling formula. Hence 

/fl K, 1/3x1 K a 

( 32 ) E (%) <^ +2<i+1+2i l^ 3 G(*)| < E (^r) ci^iCjO'CAO^ 

i=i J ' 3=1 3 ' 

for a suitable C2 > 1, where we used Y2j=[ 1 = I /?i | — 1 < 2^ 1 L 
For \/3i\ > 1, the estimate pip can then be controlled by 

(33) \(z) 2d+1 d l3l e 2ni ^^)( z )d f32 a(z + ( Vl ,u 2 ))d^G(z)\ 

< M(v 1 ,n 2 )Cl /3l l + l /32 l + l /33 l(/3 1 !/3 2 !/3 3 !) s < M(v u u 2 )C^(al) s . 

for a new constant C > 1. An easier argument shows that the same estimate holds for 
£1 = too. 

Finally, the desired result (I28p is obtained by using the estimate (I33p in (I29p . together 

with Eft+ftj+fo^ frlft-W = 3H • D 

We now show two immediate byproducts of the above theorem, namely, exponential 
sparsity of the Gabor matrix representation of T and the continuity on the Gelfand- 
Shilov spaces. Let therefore Q(g,A) be a Gabor frame for L 2 (IR d ), with g £ S^',^(M. d ), 
s > 1. Under the assumptions of the previous theorem when M = 1, we have therefore 
the estimates 

(34) |(T< 7u ,^)|<exp(-e| ? ;-x(n)| 1 / s ), 

valid for u, v G M 2( , in particular for m,!)£A. 

It is easy to see that this implies the sparsity for the Gabor matrix in the classical 
- i.e. superpolynomial - sense; cf. [U [25]. Actually, here we obtain a sparsity of 
exponential- type, as detailed in the following result. 
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Proposition 3.4. Let Q(g, A) be a Gabor frame for L 2 (R d ), with g € S s J,^(R d ), s > 1. 
Under the assumptions of Theorem \3.3l with M = 1, the Gabor matrix {Tgx^g^) is 
sparse in the following sense. Let a be any column or raw of the matrix, and let \a\ n 
be the n-largest entry of the sequence a. Then \a\ n satisfies 

Mn < Cexp (-en 
for some constants C > 0, e > 0. 

Indeed, this was shown in detail in |13|. Proposition 4.5] for any matrix satisfying an 
estimate of the type (j34]l . 

Another consequence of Theorem 13.31 and the characterization (115f) is a continuity 
result on Gelfand-Shilov spaces. 

Proposition 3.5. Let s > 1, and consider a symbol a € C 00 (M 2d ) and a phase $ 
satisfying the assumptions (fl~6|) . (fTT|) and (fTKj) with M = 1. Then the corresponding 
Fourier integral operator T in bounded on Sg(M. d ). 

Proof. The proof is analogous to the corresponding result for pseudodifferential oper- 
ators obtained in \13\ Propositions 4.7] (but here we restrict to s > 1). In short: from 
the inversion formula (|13p . we get 

V g {Tf)(v)= I (Tg u ,g v )V g f(u)du, 

JW 2d 

with g(x) = e~2< x < , say. The estimate (f34"|) together with the characterization in ^5|) 
then give the desired conclusion. □ 



4. A COUNTEREXAMPLE TO SUPER- EXPONENTIAL DECAY 

In this section we show that there is not a reasonable extension of Theorem 13.31 to 
the case s < 1. In other terms, ultra-analytic phases and symbols generally do not 
give super-exponential decay in (|24p . even for ultra-analitic windows. 

Consider, in dimension d = 1, any real- valued function (p(x), i£l, satisfying the 
following estimates: 

(35) \<P (a) (x)\ < C |Q|+1 (a!) s , Va>2, 

for some s < 1 (e.g. (p(x) = cosx). Let T be the FIO with phase &(x,n) = xr\ + (p(x) 
and symbol a = 1, therefore Tf(x) = e 271 ^^ f (x) , xiUiV) = {ViV + V(p(y)). Observe 
that the assumptions (fT6l) . (fT7|) are fulfilled, as well as (fTHl) with M = 1. Then, the 
following holds true. 

Proposition 4.1. For i/ie above operator T, suppose the following estimate holds for 
some 1/2 < s' < 1, g G 5|,'(M) \ {0}, e > 0: 

(36) <exp(-e\v- X (u)\ 1/s '), u,v£R 2 , 
Then (p(x) is a polynomial of degree at most 2. 
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Proof. The estimate (f36|) implies that if / € Sp(M.) then Tf G Sy (R) (see the proof of 
Proposition [33] or [El Propositions 4.7]). Let now f(x) = e~ x2 G S^W C Sf,'(R); 

then T/(x) = e^^e"* 2 G S^(R). The hypothesis §E§ with s < 1 and Cauchy's 
estimates imply that <p{x) extends to an entire function (p(z), z = x + iy G C. By 
Theorem 12.31 the function e 2m v( z ) e - z satisfies the growth estimate 



for some constants C, c > 0, with \x = 1/(1 — s'). The left-hand side is equal to 
e -2irim<p{z)- x 2+y^ hence 

-1m<p(z) < C(l + x 2 + |y|^). 
for a new constant C > 0. A similar estimate holds with —ip in place of (p, because — <p 
satisfies the same assumptions as ip and e - 27Tl <fi( x ) e - x g 5"|, (M) too. Therefore we get 

\1mtp(z)\ < C(l + x 2 + \yf). 

So, Im (p(z) has at most an algebraic growth, and the same must hold for the real part 
Reip(z), by the Cauchy-Riemann equations. As a consequence, \<p(z)\ has at most an 
algebraic growth, therefore (p(z) is a polynomial by the Liouville theorem. Since the 
second derivative <p"(x) is bounded by (|35j) . tp(x) must have degree at most 2. □ 

The above result shows that there is no hope to obtain super-exponential decay 
except for quadratic phases. Indeed, T is then a metaplectic operator and for those 
operators we are able to obtain optimal estimates for the corresponding Gabor matrix 
decay, as explained in the following section. 

5. A CLASS OF GENERALIZED METAPLECTIC OPERATORS 

We will study the class of Fourier integral operators whose canonical transformation 
is a linear transformation x( z ) = Az for some invertible matrix A G GL(2d, R). Since 
X must preserve the symplectic form (assumption (ii)), A must be a symplectic matrix, 
i.e. an element of the symplectic group 

Sp(d,R) = {Ae GL(2d, R) : *AJA = J} , 

where 





J ~ \h 
For z = (x,£) we shall also write 

7T( Z )f = M(F x f. 

Given A G Sp(d, R), the metaplectic operator fi(A) is defined by the intertwining 
relation 

(37) Tt(Az) = c A ^(X)7r(«)/i(^) _1 Vz G R d , 

where C4 G C, \cj\\ = 1 is a phase factor (for details, see e.g. [2 lj ) . 

If X = A = (i G Sp(d, R) , then (x, £) = (Ay+Br?, Cy+Dr?) and det || (y , t?) = 
det ^4, so that the condition ([22]) becomes det A ^ 0. 
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Viceversa, to every matrix A = G Sp(d, R) with det A ^ corresponds 

a metaplectic operator ^(A) which is a Fourier integral operator of the type (??), as 
proved in Theorem 4.51 and subsequent Remark 2 of [21], recalled below. 



Theorem 5.1. Let A = G Sp(d,R). If det A ^ we 



have 



(38) mWC*) = (det A)" 1 / 2 y e^^/fa) 

(39) $(x, 17) = -iCA _1 x + r]A~ x x - ^r]A' 1 Bi]. 

Solving ([2"U]) for the phase function in ([3U]) we obtain % = A, as expected. 

Observe that the phase in ([39]) satisfies conditions (I16D and (|17j) and the symbol 
a = (det A)" 1 / 2 in (jMJ) fulfills (pi]) with exponent s = (so also for s = 1/2) and the 
weight M = 1. So these metaplectic operators satisfy the assumptions of Theorem 13.31 
for s > 1 but the best decay result would be 

(40) \(K-A)gu,9v)\ <exp(-e|«-x(u)|), u,v £ M 2d , 

provided g G Sy^(M. d ). This decay result in not optimal, as shown by the following 
motivating example. 

Example 5.2. Consider the Cauchy problem for the harmonic oscillator: 

{.du 1 . . l2 

i— —An + 7T \x\ u = 
dt Att 11 
u(0,x) = 1*0(2:), 

with (t,x) 6Rx R rf , d > 1. For every fixed t, the solution: 

u{t, x) = (cos t)- d ' 2 [ e 2 ^^ + ^ x2+ ^f(v) dv, t ? ~ + fevr, A; G Z 

1 tan t 

can be seen as a FIO of type ([TJ with phase $t(x,r?) = zr? H (x 2 + rf) and 

cos t 2 

symbol ot = (cosi)~ rf / 2 . The associate canonical transformation is 

xt{y,v) = 



(cost)/ (— sint)7 \ f U 
(sint)7 (cost)/ J \ V 

With g{x) = e~2^ , an explicit computation shows the Gaussian decay 

\(u(t,-)g u ,g v )\ <2"f exp(-|| U - Xt (w)| 2 ), V«,t;eR 2d . 

More generally, consider the case of a FIO T with phase $ in (|39p and symbol <r that 
satisfies (fT8|) . ([19]) . that generalizes the classical metaplectic operator above, having a 
non-constant symbol. 
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Theorem 5.3. Let s > 1/2, consider a FIO T with phase $ in (|39H and symbol a 
that satisfies (I18p . (]19p . Assume g € ). T/ien i/iere exists e > suc/i t/iot 

| (Tg u ,g v )\ < M(v x , u 2 ) exp ( - e|t; - , Vu, t> E M M 

Proof. The proof uses the same pattern of the one of Theorem 13.31 The matrix of the 
second order derivatives of the phase $ is 

and the phase remainder (|26p becomes 

So |^($ 2jK)U2) (z))| < C^Hfl for 1 < |/3| < 2, whereas ^($ 2 , (vi , U2) (*)) = for 
every z £ M 2rf when |/3| > 2. In this case, by the Faa di Bruno formula, the estimate 
(|30l) is replaced by 

\Pm i \2j— Iflil a i 

IsAe^^.f^^iyW — y ft- |^| > i, 

^ ?! ^ 71 . ..7,- 

j = l J 7l + -.-+7j=/3l /J - ' J 

l<| 7fe l<2 

namely 

\dfc e 2ni * 2 >^2)M\ < cf 11 ^2 —{z) 2j ^ 11 

i=i J ' 

for a new constant C\ > 0. Now we have g € 5f so that G := g <g> g G 
and the analog of formula (I32p is here 

^^i <jr >^+M+i+^-l/»xl|5*G(z)| < ct l+l/?3| £^r(2j " WW 
i=i J ' i=i J ' 

<C7f 1+1^1(^,^1). 

where we used 

i! - j! ~° 4 Iftp 31 " 4 • 

(The last inequality holds because j < \f3i\ and s > 1/2). 
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